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Fermi-Surface Topology and Pairing Symmetry in BiS2-Based Layered Superconductors 
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In order to clarify superconducting properties of BiSi-based layered superconductors LaOi_j:Fj^BiS 2 , we evaluate the superconduct¬ 
ing transition temperature and analyze the gap function in a weak-coupling limit on the basis of the two-dimensional two-orbital 
Hubbard model. It is found that becomes maximum at x - 0.6 in the present calculations. For x < 0.45, we obtain the cZ-wave 
gap of which line nodes do not cross the pocket-like Fermi-surface curves with the centers at X and Y points. On the other hand, 
for X > 0.45, the extended s-wave gap is found for a couple of large Fermi-surface curves with the centers at M and X points. The 
variation of the Fermi-surface topology plays a key role for the gap symmetry in BiSa-based superconductors. 
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1. Introduction 


Since the discovery of new BiS 2 -based lawred superconduc¬ 
tors LaOi_jcFj,BiS 2 by Mizuguchi et aZ.yl 121 [sl], both experi¬ 
mental and theoretical investigations for the normal and super¬ 
conducting properties on BiS 2 -based materials have been per¬ 
formed intensively. This material group has the characteristic 
layered structure, similar to high-Tc cuprate superconductors, 
where denotes the superconducting transition temperature. 
Namely, BiS 2 superconducting layer is sandwiched by insulat¬ 
ing block layers. 

Among several BiS 2 -based compounds, first we pick up the 
mother compound LaOBiS 2 . This material is insulating, but by 
substituting F for O, we can dope electrons into the BiS 2 layer. 
Then, the system becomes metallic and at low temperatures, su¬ 
perconductivity occurs eventually. In particular, at x - 0.5, the 
smerconducting transition temperature Tc becomes maximum 
yt]. The highest among BiS 2 material groups has been ob¬ 
tained in the sample synthesized under high pressures |01 ■ It has 
been reported that the onset Tc is 11.1 K and the temperature at 
which the resistivity becomes zero is 8.5 K. 

Although the mechanism of superconductivity in BiS 2 -based 
materials has not been confirmed yet, the canonical model for 
this system has been proposed by Usui et al. Hi , just after the 
discovery of BiS 2 -based superconductors. The minimal model 
to describe electronic structure of BiS 2 layer is the two-band 
Hamiltonian composed of Bi 6px and 6py orbitals on the two- 
dimensional square lattice. The appearance of superconductiv¬ 
ity has been discussed from various theoretical vie^oints on 
the basis of this two-orbital Hubbard model III Isj Ip .EIH- 
The effects of electron-phonon interaction [ 13 II 14 I] and spin- 
orbit coupling II 15 I] have been also discussed. The relation be¬ 
tween the characteristic change of the Fermi-surface topology 
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and the symmetry of the superconducting gm function have 
been pointed out in previous works fiSliilllll. 


Concerning the gap symmetry, some theoretical studies sug¬ 
gested that the gap function in BiS 2 -based layered supercon¬ 
ductors could be explained by extended i-wave, tZ-wave, triplet 
or other unconventional paring scenarios ilHIES 113[H- On 


the other hand, i-wave pairin g du e to electron-phonon interac¬ 
tion has been also discussed 113 III Si. From experimental 
viewpoints, the results consistent with s- or extended i-wave 
pairing have been found in the temperature dependence of the 
superfluid density in Bi 404 S 3 fill . LaOo.5Fo.5BiS2 d], and 
NdOi-j:FiBiS 2 lull. In particular, for NdOi_^F;cBiS 2 , it has 
been reported that the gap symmetry is consistent with s- or ex¬ 
tended s-wave both for x - 0.3 and 0.5 Olflll . It is interesting 
to consider a way to reconcile those results on the symmetry of 
the gap function from a viewpoint of Fermi-surface topology 
characteristic to multi-band systems. 


In this paper, we analyze the two-dimensional two-band 
model including intra- and inter-orbital Coulomb interactions 
within a random phase approximation (RPA) to derive the ef¬ 
fective pairing interaction to induce the Cooper pair. Then, we 
solve the gap equation to obtain Tc and analyze the gap sym¬ 
metry by focusing on the Fermi-surface topology. For x - 0.3 
and 0.4, we obtain the tZ-wave gap, but the disconnected Fermi- 
surface curves with the centers at X and Y points do not cross 
the line nodes of the tZ-wave gap. When we increase x, we 
observe the elevation of except for a narrow region around 
at X - 0.45 and the conversion of the gap symmetry from d- 
wave to extended s-wave. In fact, takes the maximum value 
at X = 0.6. For x = 0.5 and 0.6, we obtain the extended s- 
wave gap for a couple of large Fermi-surface curves with the 
centers at F and M points. It is proposed that the conversion 
of the gap symmetry is characteristically induced by the varia¬ 
tion of the Fermi-surface topology in the BiS 2 -layered materi¬ 
als. Throughout this paper, we use such units as Zi = = 1. 


Preprint submitted to Journal of Magnetism and Magnetic Materials 


July 21, 2015 











2. Model and Formulation 

2.1. Model Hamiltonian 

The model Hamiltonian is given by Si 

k<TTT' ir 

+ nnna + J Yj 4<r4<r'C<-i^r'Ci2o- (1) 

i i,(T,cr' 

l.T^r' 

where Ck^icr is the annihilation operator of electron with wave- 
vector k and spin cr in the orbital t (=1 and 2), SkTT' is the 
electron energy between r and t' orbitals, n^a- - ct Qro-, 
«iT = Tio-ttiTo-, and the coupling constants U, U', J, and J' 
denote the intra-orbital Coulomb, inter-orbital Coulomb, ex¬ 
change, and pair-hopping interactions, respectively. Note the 
relations of U = U' + J + J' and J - J'. 


2.2. Electron dispersion 

As for the electron dispersion relation, by following the re¬ 
sults in the paper of Usui et aL||3l, we obtain 

e*:ii = to-I-2fi(cos^x + cos^^) 

-H 2f3 cos(kx - ky) + 2t4 cos{kx + ky) 

+ 2t(,[cos{2kx + ky) + cos{kx + 2ky)] 

+ 2ts[cos(2kx - ky) + cos(kx - 2ky)], ( 2 ) 

Sk22 - to + 2fi(cOS kx + COS ky) 

+ 2f3 COS(^;c -H ky) + 2t4 COS{kx - ky) 

+ 2t6[cos(2kx - ky) + cos{kx - 2ky)] 

+ 2ts[cos(2kx + ky) + cos(kx + 2ky)], (3) 

and 

Skl2 - £*21 

= 2f2(cos kx - COS ky) + 2f5(cos 2kx - cos 2ky) 

+ 4f7(cos 2kX cos ky - cos kx cos 2ky). (4) 

The values of to~tg are given by fo=2.811, fi=-0.167, 
f2=0.107, f3=0.880, f4=0.094, f5=-0.028, f6=0.014, tj^O.020, 
and fg =0.069 in the unit of eV Js}]. 

The band dispersion is given by 


Ek± - 


£*11 + £*22 ± -^(£*11 - £* 22 )^ + £] 


.2 

*12 


(5) 


Hereafter we define Ek+ = E^a and = E^b, which are 
called the A and B bands, respectively. 

The band structure has been discussed in Ref. j^l, but in order 
to make this paper self-contained, we show E^a and E^b along 
the path with high symmetry on the Brillouin zone in Fig. 1(a). 
For X - 0.3 and 0.4, as shown in the insets, we observe that 
the Fermi levels cross only the B band near the X point, while 
between M and F points, they do not cross the band B. For 
X = 0.5, the Fermi level turns to cross the band B between M 





Figure 1: (a) EkA and Eks along the path of F—>X^M—»r. The horizontal 

broken lines denote the positions of Fermi levels for x = 0.3, 0.4, 0.5 and 0.6 
from the bottom to the top. Insets denote the band structure in magnified scales 
around at (tt, 0) and a point between (tt, n) and (0,0). Fermi-surface curves for 
(b) X = 0.3, (c) X = 0.4, (d) x = 0.5, and (e) x = 0.6. 


and F points, but it is still under the band A at the X point. For 
X = 0.6, the Fermi level crosses not only the B band but also 
the A band near the X point. 

The above characters appear in the change of the Fermi- 
surface structure, as shown in Figs. 1(b)-1(e). In Fig. 1(b), we 
depict the Fermi-surface curves for x = 0.3. A characteris¬ 
tic issue is that pocket-like disconnected Fermi-surface curves 
appear around at the X and Y points. We remark that these 
Fermi-surface curves originate from the B band and they do not 
cross the lines of kx - ±ky. This issue will be important in 
the discussion of the gap function later. Another characteristic 
point is the shape of the Fermi-surface curves. We observe the 
disconnected Fermi-surface curves, which possess the projec¬ 
tions extending to the direction to those of other Fermi-surface 
curves. For x - 0.4, we still observe pocket-like disconnected 
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Fermi-surface curves around at the X and Y points, as shown in 
Fig. 1(c). Note that the projections further extend to the direc¬ 
tions to other Fermi-surface curves, but they do not touch. 

When we further increase x, the projections eventually touch 
other Fermi-surface curves and the disconnected structure is 
dissolved, as observed in Fig. 1(d) for x = 0.5. Then, we ob¬ 
serve a couple of large Fermi-surface curves with the centers at 
the F and M points. These two Fermi-surface curves originating 
from the B band are formed by the merge of the projections of 
pocket-like Fermi-surface curves in Fig. 1(c). We should pay at¬ 
tention to the existence of the Fermi-surface curves on the lines 
of kx - ±ky. We also remark that the topological change of the 
Fermi-surface structure occurs around at x = 0.45, as already 
mentioned in Refs. ll3l and M- 

Now we further increase the values of x. The Fermi-surface 
curves for x = 0.6 are shown in Fig. 1(e). Flere we observe 
the large Fermi-surface curves around at the F and M points, al¬ 
though the shape is slightly deformed from those for x = 0.5 
in Fig. 1(d). Note also that new pocket-like Fermi-surface 
curves appear around at the X and Y points. These curves orig¬ 
inate from the A band, in sharp contrast to other Fermi-surface 
curves. We find that there occurs the topological change for the 
appearance of the Fermi-surface curve originating from the A 
band at x = 0.52, as already pointed out in Ref. lliol] . 


2.3. Spin and orbital susceptibilities 


In order to investigate superconductivity yound the spi n 
and/or orbital ordered phases 2j., 22, 3 ^ 26,271. 


we calculate spin and orbital susceptibilities, x^{q) and x°(q), 
respectively. Here we follow the formulation and the notations 
in Ref. yOfl. Within the RPA, spin and orbital susceptibilities 
are giyen in the 4x4 matrix form as 


r(q) - [1 - ( 6 ) 


and 

riq) = [i + U°x(q)r'x{q), (7) 


respectiyely. Here labels of row and column in the matrix ap¬ 
pear in the order of 11, 22, 12, and 21, which are pairs of orbital 
indices 1 and 2. Note that 1 is the 4x4 unit matrix. The interac¬ 
tion matrices W and U° are giyen by 




U J 
J U 
0 0 
, 0 0 


0 

0 

U' 

r 


0 'I 
0 

r ’ 

U' 


( 8 ) 


and 




U 2U' -J 
2U'-J U 

0 0 

0 0 


0 0 1 

0 0 

2j - u' r 
r 2J-U' , 


(9) 


respectiyely. Each matrix element of xiq) is defined by 

X,v.ali{q) = -7’ 2 iUn), (10) 

k,n 


where T is a temperature, G^^Xk,ia>„) is the non-interacting 
Green’s function propagating between p- and v-orbitals, and 
Un - nT{2n H- 1) with an integer n. The instabilities for the 
spin- and orbital-ordered phases are determined by the condi¬ 
tions of det[i - f/’*;f'(^)]=0 and det[i + lJ°x{q)\=0, respectiyely. 

2.4. Gap equation 

By using;f'®(^) wdx°{q), we obtain the superconducting gap 
equation in a weak-coupling limit as 

A^(jt) = ^ vHk - (11) 

k' 

where A.^{k) =[Ajj(k), A 22 (k), A 2 j(k)]' is the gap func¬ 

tion in the yector representation for singlet (^ = S) or triplet 
(^ = T) pairing state. The matrix elements of the singlet- and 
triplet-pairing potentials are, respectiyely, giyen by 

( 12 ) 

(13) 

which are composed of spin and orbital susceptibilities as 

W\q) = 0 * + LiY{q)u\ (14) 

W°(q) = -tj° + u°riq)u°. (15) 

The element of the pair correlation function ^{k) is giyen by 

= r ^ GW(k, -iaj„)Gf^(-k, ioJn). (16) 

n 

We obtain by solying the gap equation eq. (fTTT i. Note that 
Tc is defined as a temperature at which the positiye maximum 
eigenyalue of the gap equation becomes unity. For the eyalua- 
tion of the eigenyalue of eq. (fTTb . we use the power method. 
The first Brillouin zone is diyided into 64 x 64 meshes and 
we exploit the fast-Fourier-transformation algorithm, in order 
to accelerate the sum of large amount of momenta. 

As for the yalues of Coulomb interactions, we note the rela¬ 
tions of t/ = U' + J + J' and J - J', as mentioned aboye. Here 
we assume the yalue of 7 as 7 = t//6. Then, from the aboye 
relations, we obtain U' - 2t//3. In this paper, we fix the yalue 
of 1/ as t/ -2 eV. Note that the energy unit is set as eV in 
this paper. The bandwidth W of the present tight-binding band 
structure is W = 4.5 eV. Thus, the yalue of 1/ = 2 is not so large 
in comparison with W. 

3. Calculation Results 

3.1. Phase diagram 

Let us show our calculation results for the transition tem¬ 
perature Tc as a function of doping x in Fig. 2. First of all, 
we should note the magnitude of the yertical temperature axis. 
Since the energy unit is eV in the present calculations, the max¬ 
imum yalue of becomes about 600 K, if we simply use the 
present energy unit. Here we should pay attention to the fact 
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Figure 2: Phase diagram on the plane of doping x and temperature T. Solid 
black circles denote the magnetic instability temperature I’m, while solid red 
squares and blue diamonds indicate Tc for the extended s- and t/-wave symme¬ 
tries, respectively. Note that the solid curve denotes the phase boundary, but 
the broken curve does not. Inset shows the part of the phase diagram around at 
X = 0.45 in a magnified scale. 

that we have not included at all the effects of the retardation of 
effective interactions and the quasi-particle life-time. If we cor¬ 
rectly include the vertex corrections in addition to above strong¬ 
coupling effects, Tc is considered to be strongly suppressed. 
Thus, it is not necessary to consider seriously the magnitude 
of Tc in the present calculations. Rather we should pay our at¬ 
tention to the relative change of Tc by the doping x. 

The lower curve denotes the magnetic instability temperature 
I'm determined from det[l-0^®;f'(^)] = 0. In the present calcula¬ 
tions, when we decrease the temperature, the spin susceptibility 
always diverges prior to the orbital one and we do not consider 
the orbital susceptibility. Then, we solve the gap equation in 
the range of T > Tm to consider the superconductivity induced 
by magnetic fluctuations in the vicinity of the magnetic insta¬ 
bility temperature. Note that we cannot discuss the coexistence 
of magnetism and superconductivity in the present RPA. 

In Fig. 2, we immediately notice that the superconducting 
phase in the weak-coupling limit appears in the wide range of 
the values of doping x and temperature T. Roughly speaking, 
the shape of the curve for is similar to that for I'm. In par¬ 
ticular, the position of the maximum is exactly the same as 
that of I'm. We also point out that the superconducting phase is 
divided into two regions: The d-wawe phase for x < 0.45 and 
the extended i-wave phase for x > 0.45. The figure around at 
X 0.45 in a magnified scale is shown in the inset of Fig. 2. 
In the present calculations, in a very narrow region around at 
X = 0.45, Tc is suppressed. Instead, we find the ferromagnetic 
(FM) phase in this region, since the ordering vector is found 
to be ^ = (0,0) in this region. Near x = 1.0, Tm becomes 
zero since both U and the density of states at the Fermi level 
are small in this case, but Tc keeps a large value, although it is 
smaller than the peak value at x = 0.6. For x < 0.2, we cannot 
obtain both magnetic and superconducting phases in the RPA 
for the system with small electron densities. 


We note that there is not significant structure in around 
at X = 0.52, at which the Fermi-surface topology is changed 
due to the appearance of the Fermi-surface curves originating 
from the A band. We deduce that the effect of the pairing in the 
A band will not be so significant, since the superconductivity 
occurs mainly in the B band. In the d-wawe phase, except for 
the narrow region around at x = 0.45, Tc is increased with the 
increase of x. On the other hand, in the extended s-wave phase, 
Tc takes a broad peak at x = 0.6. Note that the magnetic in¬ 
stability curve also exhibits the peak at x = 0.6, but the peak 
structure is rather sharp. 

Let us now compare our phase diagram with experimental 
results. First we briefly review the experimental facts. For 
LaOi_xFjcBiS 2 , the doping dependence of Tc seems to be sen¬ 
sitive to the way of sample synthesis For as-grown 

samples, Tc is increased between x = 0.2 and 0.5, whereas it 
becomes almost constant between x = 0.5 and 0.7. For an¬ 
nealed samples under high pressures, Tc is totally increased and 
it becomes maximum at x = 0.5. In addition, Tc seems to be 
abruptly increased between x = 0.4 and 0.5. 

For NdOi^xFjBiS 2 , we have found weak doping dependence 
ofTc 0. It has been found that superconductivity appears in 
the region of 0.1 < x < 0.7 and Tc becomes maximum at x = 
0.4, although the change of Tc by doping x is not so significant 
in comparison with LaOi_xFxBiS 2 . 

In our phase diagram, the maximum Tc is found at x = 0.6, 
which is deviated from x = 0.5 for LaOi_,iFj^BiS 2 and x = 
0.4 for NdOi_jcFj:BiS 2 . This deviation may be explained by 
the strong coupling calculation or by the further inclusion of 
electron-phonon interaction. These points will be discussed in 
future. We remark a point that the nominal value of x may be 
deviated from the amount of electron doping in BiS 2 layer 113 ill . 
We should pay due attention to this point in the comparison 
with experiments. It is another future task. 

Around at x = 0.45, the Fermi-surface topology is changed 
and the singlet superconductivity has been found to be sup¬ 
pressed. As expected from the FM phase around at x = 0.45, 
the eigenvalue of the triplet pair has been found to be larger 
than that of the singlet one in this region, although it has not 
reached the unity in the present calculation. The tendency of the 
triplet pair may be also related to the previous work in Ref. . 
When we include the strong-coupling effects and the vertex cor¬ 
rections in the evaluation of the normal and anomalous self¬ 
energies, the triplet superconductivity may occur around at 
X = 0.45. This is one of interesting future issues. 

3.2. Susceptibilities and Fenni-surface curves 

In order to grasp the kind of the magnetic fluctuations for the 
pair formation, we show the spin susceptibilities in the RPA for 
X = 0.3, 0.4, 0.5, and 0.6 in Figs. 3 and 4 along the path of 
q - (0,0)^(7r, 0) —»( 2 r, 7r)^(0,0). Since we have found that 
orbital susceptibilities are totally suppressed and they do not 
show any significant structures in comparison with spin ones in 
the present calculations, we show only spin susceptibilities in 
this paper. 

In Figs. 3(a) and 3(b), we show the spin susceptibilities, re¬ 
spectively, for X = 0.3 and 0.4. Note that the temperatures are 
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Figure 3: Spin susceptibilities in the RPA for (a) x = 0.3 and (b) x = 0.4. Inset 
in each panel denotes the Fermi-surface curve. 


Figure 4: Spin susceptibilities in the RPA for (a) x = 0.5 and (b) x = 0.6. Inset 
in each panel denotes the Fermi-surface curve. 


set as Tc, which are - 0.0255 for x — 0.3 and - 0.0355 
for X = 0.4. We hnd the enhanced spin fluctuations of q=(q, q), 
which will contribute to the Cooper-pair formation for super¬ 
conductivity. The value of q becomes small when x is increased 
from 0.3 to 0.4. For x - 0.4, we observe another small peak for 
large q. Note that the FM fluctuations also grow for x - 0.4 
in comparison with that for x - 0.3. The FM fluctuations tend 
to be significant, when x approaches 0.45, as expected from the 
existence of the FM phase around at x = 0.45. 

In Fig. 4(a) and 4(b), we exhibit the spin susceptibilities in 
the RPA for x - 0.5 and 0.6, respectively. Note that the tem¬ 
perature is set as T = 0.0445 and 0.0565, which are equal to 
at X = 0.5 and 0.6, respectively. In the bare susceptibility, the 
maximum peak appears at ^ = {q, 0) with small q. This peak 
structure is due to the nesting of a couple of large Fermi-surface 
curves with the centers at the F and M points fl. In fact, for 
X = 0.5, the peak at ^ = {q, 0) is larger than other peak values in 
the RPA spin susceptibility, although anti-ferromagnetic (AF) 
spin fluctuations also grow. For x = 0.6, the peak at ^ = {q, 0) 
is still signihcant, but we observe a clear peak at ^ = (tt, n), 
suggesting AF spin fluctuations. We deduce that the evolution 
of the AF spin fluctuations is closely related to the elevation of 
Tc from X = 0.5 to x = 0.6. 


3.3. Gap functions and Fermi-surface curves 

Now we show the results of the gap function for the Cooper 
pair from the same band, Apfk) and Ab(^), which are ob¬ 
tained by the appropriate linear combinations of Aj j (k), 

Aj 2 (k), and Ajj)^). There exist the components of the Cooper 
pairs formed by different band electrons, but we do not show 
them here, since their magnitudes are small in comparison with 
those of the same band. Note also that only the singlet pair has 
been stabilized, although we have investigated both singlet and 
triplet channels. Thus, we show here the results for singlet gap 
functions. 

In Figs. 5, we exhibit Apf k) and Ab)^) for x = 0.3 and 
T - 0.0255 by using the color gradation to visualize the 
momentum dependence at a glance. First we notice that the 
gap functions Apfk) and Ab(^) possess the ^f-wave symmetry, 
which is characterized by the sign change of the gap for the in¬ 
terchange of kx and ky. In this doping, there is no Fermi-surface 
curve from the band A, but we observe the signihcant ampli¬ 
tude around at the X and Y points, since the band A is near the 
Fermi level on these points, as shown in Fig. 1(a). 

As for the gap function on the band B in Fig. 5(b), we ob¬ 
serve that the nodes of the d-wave, kx - ±ky, are running be¬ 
tween the pocket-like Fermi-surface curves. Namely, the nodes 
do not cross the Fermi-surface curves. However, we hnd the 
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Figure 5: (a) Aa(*:) and (b) Ab(A:) for = 0.3. We show the Fermi-surface Figure 6: (a) Aa(A) and (b) Ab(*) for x = 0.4. We show the Fermi-surface 
curves in the figure. curves in the figure. 


sign change of the gap at the heads of the projections for the 
interchange of and ky. Note that the signs of the gap func¬ 
tions near the X and Y points are also different from those at 
the heads of the projections of the Fermi-surface curves. For 
X - 0.3, we can observe other node lines connecting the points 
of {-7T, 0), (0, -7t), {tt, 0), and (0, n). These node structures seem 
to be consistent with those proposed by Usui et al. ||3l. 

In Figs. 6, we show Aa(^) and Ab(^) for x - 0.4 and 
T - 0.0355. The tf-wave node structure is similar to that in 
Figs. 5, but we do not observe the node lines connecting four 
points of (-TT, 0), (0, -n), (tt, 0), and (0,7r) in Fig. 6(b). Rather 
the gap function Ab(^) changes its sign for three times on the 
straight line from (tt, 0) to (0, tt). The gap function Ab(^) seems 
to includes the components of the symmetry higher than the 
t/-wave, for instance, g-wave Eh. 

The gap structure on the Fermi-surface curves around at the 


X and Y points has been investigated in NdOo. 7 iFo. 29 BiS 2 with 
the use of an angle-resolved photoemission spectroscopy lli^l . 
The large anisotropy has been observed and the gap size on the 
point at which the Fermi-surface curve crosses the line from 
(tt, 0 ) to (tt, tt) is smaller than that on the point at which the 
Fermi-surface curve crosses the line from (tt, 0) to (0,7r). The 
tendency seems to agree with the present results, but in order to 
conclude the existence of the node, it is necessary to perform 
further experimental and theoretical investigations. 

Next we consider the gap functions for x > 0.45. In Fig. 7(a) 
and 8(a), we depict the gap functions A^ik) for x = 0.5 and 0.6, 
respectively. Note that for x = 0.6, we also depict the Fermi- 
surface curve originating from the A band around at the X and 
Y points. In both cases, we do not find the nodes characteristic 
to the d-wa\e in the gap function. The node lines satisfying 
cos kx + cos ky = 0 seem to be observed, but it is difficult to 
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Figure 7: (a) AaC*) and (b) Ab(A:) for jc = 0.5. We show the Fermi-surface Figure 8: (a) Aa(A) and (b) Ab(*) for x = 0.6. We show the Fermi-surface 
curves in the figure. curves in the figure. 


confirm them only from this hgure. 

In Figs. 7(b) and 8(b), we show the gap functions Ab(^) and 
the large Fermi-surface curves for x = 0.5 and 0.6, respectively. 
Also in these cases, the nodes characteristic to the li-wave gap 
cannot be observed. Note that the structures in the gap functions 
are quite similar to each other irrespective of the difference in 
X. We hnd several numbers of node lines expressed as ky = 
±kx + C, where C indicates an appropriate constant. When we 
change k in Ab(^) along the line from the F to M points, the 
gap changes its sign for several times. 

We remark that the Fermi-surface structure becomes differ¬ 
ent between x - 0.5 and 0.6. In particular, the curvature of the 
Fermi-surface curves for x - 0.5 is larger than that for x - 0.6. 
Since the nodes are well expressed by the lines of ky - ±kx -i- C, 
the nodes of the gap easily appear on the Fermi-surface curves 
with large curvature for x = 0.5 in comparison with those for 
X = 0.6. From Figs. 7 and 8, there seems to occur a peculiar 


possibility that the Fermi-surface curves tend to avoid the ex¬ 
tended i-wave gap functions with large amplitude. This issue 
should be investigated in detail in future. 

4. Discussion and Summary 

In this paper, we have solved the gap equation to obtain the 
superconducting transition temperature and the gap function 
within the RPA on the basis of the two-band Hubbard model 
for BiS 2 -based layered compounds. When the topology of the 
Fermi-surface curves of the B band drastically changes around 
at X = 0.45, we have observed the change of the symmetry 
of the gap function. For x < 0.45, we have obtained the d- 
wave gap with the nodes on the lines of ky = ±ky, which do 
not cross the Fermi-surface curves with the pocket-like discon¬ 
nected shapes around at the X and Y points. On the other hand. 
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for X > 0.45, we have found the extended s-wave gap. An inter¬ 
esting issue is that the amplitude of the gap seems to be small 
or zero on the Fermi-surface curves. 

The effect of the nodes on the temperature dependence of 
physical quantities is quite interesting for the comparison with 
experimental result. In particular, it has been reported that 
the temperature dependence of the penetration depth indicates 
the s- or extended s-wave gap both for x - 0.3 and 0.5 in 
NdOi^j:Fj^BiS 2 ifisl] . In the present results, we have found the 
t/-wave gap for x < 0.45, but the main nodes of t/-wave gap 
do not cross the Fermi-surface curve. For x > 0.45, we have 
obtained the extended s-wave gap. However, for both gap func¬ 
tions, we have found the node structure on the Fermi-surface 
curves depending on the doping x. It is interesting to clarify the 
temperature dependence of the penetration depth for the gap 
function with the present node structure. This is one of future 
problems. 

In the present calculations, we have considered only the 
Coulomb interactions. However, around at the Lifshitz tran¬ 
sition point at which the Fermi-surface topology is abruptly 
changed, the coupling between the charge fluctuations and lat¬ 
tice vibrations seems to be important. In fact, the role of 
pho nons in BiS 2 -based materials has been discussed intensively 
1 12,12 14 1. We believe that it is necessary to include electron- 
phonon interaction to the present two-orbital Hubbard model. 
It is another future problem. 

In summary, we have discussed the gap symmetry of BiS 2 - 
based layered superconductors by analyzing the two-orbital 
Hubbard mode within the RPA. For x < 0.45, we have found 
the d-wave gap on the pocket-like disconnected Fermi-surface 
curves, whereas for x > 0.45, the extended i-wave gap has been 
found on a couple of large Fermi-surface curves with the cen¬ 
ters at the F and M points. The symmetry change in the gap 
function is believed to be induced by the topological change in 
the Fermi-surface structure. 
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